DNA electrophoresis in designed channels 
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We present a simple description on the electrophoretic dynamics of polyelectrolytes going through 
designed channels with narrow constrictions of slit geometry. By analyzing rheological behaviours of 
the stuck chain, which is coupled to the effect of solvent flow, three critical electric fields (permeation 
field £ (per) ~ AT 1 , deformation field E {deS) ~ N' A/b and injection field E {inj) ~ TV , with TV 
polymerization index) are clarified. Between E^ perS> and E^ m ^ , the chain migration is dictated by the 
driven activation process. In particular, at E > E^ de ^ , the stuck chain at the slit entrance is strongly 
deformed, which enhances the rate of the permeation. From these observations, electrophoretic 
mobility at a given electric field is deduced, which shows non-monotonic dependence on TV. For long 
enough chains, mobility increases with TV, in good agreement with experiments. An abrupt change 
in the electrophoretic flow at a threshold electric field is formally regarded as a nonequilibrium phase 
transition. 

PACS numbers: 83.50.-v, 36.20.Ey, 87.14.Gg 



I. INTRODUCTION 

The fractionation of DNA molecules depending on 
their length is a required step in many sectors of bio- 
logical sciences. A standard technique for this purpose is 
a gel electrophoresis, in which negatively charged DNA 
molecules migrate driven by electric field inside the ran- 
dom environment of a gel 0, 12 ■ During the driven diffu- 
sion, DNAs have to negotiate with the geometrical con- 
straints imposed by the gel, resulting in lower mobility 
for longer chains. It is well known, however, that there 
is an upper limit for the length, beyond which mobility 
of DNAs become length- independent. 

Many efforts have been carried out to improve on this 
restriction. Among them, artificially designed devices 
with micro-patterns have been getting a lot of attention 
[J • While in gels the environment for the tar- 
get molecule is random and there is an upper bound for 
the statistical "mesh" size set by the mechanical weak- 
ness of dilute gels (for example, the maximum mesh size 
in agarose gels is ~ 0.3/zm), the geometry of the device 
with much larger size constriction can be designed in a 
desired way. Another advantage is that the microfabri- 
cated device is easy to build in integrative separation sys- 
tems, which would be required in the field of biomolcculcs 
analysis in the future. 

Cleary, the progress in this area relies on the under- 
standing of dynamics of polymers in such restrictive en- 
vironments. While the dynamics without or with weak 
external bias has been extensively studied^. ITol ITTl Il2| . 
little is known about the dynamics involving out of equi- 
librium, large deformation caused by the external field. 
Indeed, long polymers are so soft that they could be eas- 
ily deformed under weak external forces. Examples are 
a chain under stretching 13] and a tetherd chain under 
hydrodynamic flowplj. etc.. The same also applies for 
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the polymer dynamics in restrictive environments, the 
consequences of which seem not to have been explored 
well. 

In this article, we investigate the electrophoretic dy- 
namics of a polyelectrolyte chain in an artificial chan- 
nel, where two domains of slit constituting of different 
depths are alternately arranged (Fig. [JJ, with dc electric 
field. The choice of this particular geometry is moti- 
vated by recent experimental studies, which suggest the 
potentiality for the DNA separation over the wide range 
of chain length 0, [ij ■ A remarkable feature in this sys- 
tem is that longer DNAs migrate faster in contrast to 
usual situations in gel electrophoresis, which triggered 
further ex peri mental and several computer simulation 
studies [lM VWL \VA Il8| . Here we present a simple scaling 
argument on the electrophoretic dynamics, which takes 
the chain rheological behaviours and hydrodynamics in- 
teractions into account simultaneously, in the hope of 
providing lucid descriptions. An abrupt change of the 
dynamical state at a certain threshold electric field, as 
was observed in experiments, is interpreted as a nonequi- 
librium phase transition by taking an appropriate limit. 
In addition, we briefly discuss the dynamical diagram ob- 
served in the gel electrophoresis, in particular the tran- 
sition from entropic trapping to reptation regimes upon 
increasing the electric field. 

In Sec. IHJ we briefly remind of the basis of the dy- 
namics of polyelectrolytes under the electric field. Then, 
in Sec. IIIII and IIVI we analyze the permeation process 
of polyelectrolytes into the narrow channel. In Sec. E] 
we discuss the electrophoretic mobility, in particular, its 
dependence on the chain length and the electric field. We 
predict the four distinct regimes for the electrophoretic 
dynamics; with the increase in the electric field, one en- 
counters in order, (i) unbiased or weakly biased dynam- 
ics, (ii) biased dynamics, (hi) biased dynamics with the 
significant deformation of the stuck chain, and (iv) free 
migration without the trapping. If one relies on the "lo- 
cal friction" picture (neglect of the hydrodynamics), the 
regime II disappears, which is discussed in Appendix lAl 
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To highlight the effect of confinement geometries, we dis- 
cuss in Appendix [B] the electrophoretic dynamics in a 
cylindrical geometry. We also address the connection to 
the flow-injection problem, in which polymers are driven 
by the solvent flow instead of the electric field (Appendix 
0. 
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FIG. 1: A schematic picture of the devise. It consists of a 
series of wide and narrow slits, whose depth (length) are D 
(L D ), D (Ld), respectively. 



II. DYNAMICS OF POLYELECTROLYTES 
WITH GEOMETRICAL TRAPS 

Consider a polyelectrolyte chain with N segments of 
length a, each of which carrying electric charge of qe 
(e: elementary charge). For strongly charged polyelec- 
trolytes such as DNA, this charge should be regarded as 
the renormalized one due to the counterion condensation. 
As usual for electrophoresis buffers, the ionic strength 
of the solution is rather high such that the bulk chain 
size is given by that of a neutral chain in good solvent, 
R ~ aN 3 / 5 . It also affects the dynamics through the 
screening of hydrodynamic coupling on scales larger than 
Debye length, which leads to the so-called local friction 
picture, i.e., Rouse dynamics, for the free electrophore- 
sis. As a consequence, polyelectrolytes in the free so- 
lution are treated as free draining coils and the elec- 
trophoretic mobility fi e i — qe/rja (77 being the solvent vis- 
cosity) is neither dependent on the chain length nor the 
conformational]. However, as was pointed out in ref.Q, 
for polyelectrolytes interacting with an immobile obstacle 
(more generally, for polyelectrolytes under the simultane- 
ous action of electric fields and non-electric forces), this 
local friction picture breaks down. In such situations, one 
has to deal with the electrohydrodynamics, which is cou- 
pled with the deformation behaviours of the chain. Al- 
though this does not allow, in general, a simple analytic 
treatment, it is possible to linearize the coupled electro- 
hydrodynamic equations in the case relevant to typical 
electrophoresis experiments (low Raynolds number situ- 
ations). This scheme allows the superposition of electric 
and non-electric forces, which leads to a hydrodynamic- 
electric equivalence 0, 0, la- It says that, for example, 
the deformation of an end-anchored polyelectrolyte un- 
der the uniform electric field E is equivalent to that of 
an end-anchored polymer in a uniform solvent flow at 
speed v = fi e iE. And the force to immobilize that poly- 
electrolyte (stall force) f stall = l^eiE, with 7 being the 



friction coefficient for the chain in its steady-state con- 
formation. In the situation of our interest, the polymer 
may be blocked at the entrance of the narrow slit. It may 
be deformed there, but does not show any net motion, 
while the counterions are still driven without any obsta- 
cle, which induces the solvent flow. Taking these points 
in mind, in what follows, we consider the dynamics of 
polyelectrolytes in designed channels by focusing on the 
elementary step of the channel device, which consists of 
one narrow slit facing large enough space (Fig. 



III. UNBIASED DYNAMICS 

Wc first consider the unbiased or weakly biased dynam- 
ics, where the effect of electric field is very small or even 
negligible. The dynamics of polymers in restrictive envi- 
ronments depends on the relative size ratio between the 
polymer and the typical size of the confinement, which 
may be called "pore" . The situation of our present in- 
terest corresponds to the pore size comparable to, or a 
little smaller than the polymer size in a bulk. Then, 
it is known that the free motion (i.e., without external 
force) of polymers is described by the so-called entropic 
trapping process IS ES IH 03 The physics involved 
in this process is essentially the same as size-exclusion 
chromatography, which states that there is a free energy 
cost AFo(N) for the polymer to be squeezed in a narrow 
space. It should be noted that even in the situation with 
non-adsorbing wall in good solvent, the origin of con- 
finement free energy AFq(N) comes from the segmental 
interactions or the reduction in the conformational en- 
tropy depending on the architecture of the polymer and 
the geometry of the confinement, which is reflected in the 
scaling form of AF 0. Since the case of our interest, 
i.e., a linear polymer confined in a slit, corresponds to the 
entropic confinement, the terminology of entropic barrier 
would be appropriate. In this case, the scaling form of 
the confinement free energy is [IS l20| 



AFq(JV) 
k B T 



N 



(1) 



This free energy barrier leads to the exponential slowing 
down of the chain mobility with the increase in chain 
length and/or the decrease in the pore size: 



AF(N)\ 



(2) 



IV. STRONGLY BIASED DYNAMICS 

If the motion of the chain is biased by the external 
force, one has to modify the above discussion for the un- 
biased dynamics. There are several attempts to this di- 
rection including the injection of linear chains pH and 
branched chains |22| into a narrow capillary by hydrody- 
namic flow, biased electrophoretic mobility of DNA using 
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a slit channel [H, H, EH EH E3, E3- In these studies, a 
common strategy is that one notices events occurring in- 
side the channel only, which we call "inside approach" . 
The validity of the inside approach can be justified in 
some cases 0, • However, the events occurring out- 
side of the channel become important depending on the 
dissipation mechanism and/or the space dimension of the 
channel ( Appendix lAl lB|l . Indeed, we will show that the 
system under consideration falls upon the latter case. For 
the sake of convenience, nevertheless, we first present the 
inside approach. The modification due to the outside 
events is necessary at high electric field, which will be 
done later. 



(a) 






Xd 








of blobs in a slit, the primarily injected blob will proceed 
straight toward the direction of the electric field, which 
indeed corresponds to the most energy-saving strategy. 
In other words, the nucleated hernia would be more or 
less straight. This means that the essential part of the 
permeation process is described as the one-dimensional 
problem (Fig. This notion is consistent with the 

video images of ref. |5j- 

Let us denote the length of a sequence of blobs already 
permeated as y. The free energy in the corresponding 
state is 



F in (y) 



(ILD 2 - qeEi 



dy 



(3) 



The first term in the integrand comes from the confine- 
ment effect, where II ~ fcsT/(Z? 3 ) is the osmotic pres- 
sure. The second term represents the force imposed by 
the electric field inside the slit E in . Equation © indi- 
cates the presence of a free energy barrier at y = y* n : 



AR, 



!,- u T a< E;„ ^ D 2 ) 



and 



y, 



k B T /a 
qeE in \D 



(4) 



(5) 



The main result of the inside approach with the external 
force is that the free energy barrier does not depend on 
the chain length N, which is in a remarkable contrast to 
eq. (HJ. 



FIG. 2: A schematic picture of the elementary unit of the 
channel device ((a)side and (b)top views) and a process of 
polyelectrolyte permeation by electric field, where the chain 
section outside the narrow channel is assumed to be unper- 
turbed ("inside approach"). In (b) and FIG. [31 the shaded 
part designates the narrow region. 



A. Inside Approach 

We now consider the situation, where a large poly- 
electrolyte chain is facing the entrance of the narrow slit 
(Fig. |2J). Then the process of chain permeation may be 
induced by the electric field. In the early stage of the 
permeation, the insertion of a single blob of size D is re- 
alized easily. The corresponding number of monomers is 

3 /5 

denoted as go'- D = ag^ . 

But for the chain to be further permeated, a sufficient 
electrostatic force is required to overcome the osmotic 
force due to the confinement. If the moderate range of 
the electric field is applied, this permeation process is 
formulated as a barrier crossing problem Although 
there is a two-dimensional freedom for the arrangement 



B. Modification at High Electric Field: Single 
Chain Rheology 

The main assumption of the inside approach is the ne- 
glect of the chain portion still outside of the narrow chan- 
nel. However, this assumption is not always valid in the 
case of our problem, since as we shall show below, the 
chain facing the entropic barrier, i.e., the chain stuck at 
the entrance of a narrow region, is deformed to a notable 
degree under the influence of the rather high external 
field. Then, one has to deal with the single chain rheol- 
ogy; how the chain is deformed under the influence of the 
external force. This case with high external force is in- 
deed interesting in the context of the application, which 
requires the shorter separation time. 

Let us first clarify when the deformation of the stuck 
chain becomes evident. In experiments, the dielectric 
constant of the channel is usually much lower than that 
of the water solution. Then, the electric flux lines inside 
the channel interior scarcely leak out. From the conser- 
vation of electric flux and the geometrical symmetry of 
the channel, the electric field at distance r from the en- 
trance is evaluated as E(r) ~ Ei n D/r. In the light of a 
hydrodynamic-electric equivalence, we calculate the stall 
force for the stuck chain. Under the action of the electric 
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field, the chain is pressed against the wall, and its shape 
is distorted to non-spherical form. If the electric field is 
small enough, it would be plausible to assume that the 
average chain shape is non-spherical (reflecting the sym- 
metry of the electric field), but the overall size is still 
given by R ~ aN 3 / 5 (unpertubed size). Assuming that 
the stuck chain is not essentially unperturbed (overall 
chain size is given by R = aiV 3 / 5 ), the electrostatic force 
acting on the whole chain is evaluated as 



qeE in D 



(6) 



where 70 ~ r/R is the friction coefficient for an unper- 
turbed chain, and E ~ ^ E{r)r/ R 2 dr is the average 
electric field acting on the stuck chain. If this force 
is smaller than k B T/R, the assumption of unperturbed 
chain conformation is valid. In the opposite case, how- 
ever, the chain is significantly deformed due to the con- 
vergent electric field. Equating k B T/R with eq. iJBJ, we 
have the chain length dependent threshold electric field 
(in a narrow slit): 



£(<fc/) 



N 5 

qeD 



(7) 



Equation Q indicates that the longer chain is more easily 
deformed. 




Z S >D 



FIG. 3: A stuck chain at the entrance of the narrow channel 
under strong electric field (E {def) < E < E iinj) ). A chain is 
strongly deformed by the convergent field, and its conforma- 
tion is represented as a parallel alignment of blobs of size £. 
(top view corresponding to FIG. |U(b)) 



At Ein > E^ de ^ , the chain is intensely pressed against 
the slit entrance, and the molecular shape becomes uni- 
axial something like a cigar. The resulting deformed con- 
formation is described as a parallel alignment of blobs of 
size f along the slit entry (Fig. |3J). Inside blobs, the force 
due to electric field is not significant, thus, the blob size £ 
is related to the number of monomers inside it g as in the 
bulk, i.e., £ = ag 3 ^ 5 . Accordingly, bios size £ is obtained 
by the following equation 



k B T 

h 



(8) 



where is the force acting on each blob; 

f s ~ qeE ln D/a (9) 
From these, we obtain 

kftT / a 



qeE Vt 



(5) 



(10) 



Therefore, the length of the stuck chain along the slit 
entrance is 



(11) 



If the blob size £ is smaller than the slit size D, there is no 
entropic barrier arising from the difference in the osmotic 
pressure; the chain will be injected smoothly inside the 
slit. The critical field for the injection is, thus, obtained 
from eq. (|lf)|> with £ = D as 



k B T / a \ 



qe 



\D 2 ) 



D 



(12) 



At electric field above £?( OT J), the chain trapping at the 
slit entrance is irrelevant, therefore, the local friction as- 
sumption becomes valid for the description of the clcc- 
trophoretic dynamics of polyelectrolytes. Then the elec- 
trophoretic mobility is independent on both the chain 
length and the electric field. Note that the same phe- 
nomenon for the polyelectrolyte without hydrodynamic 
interactions corresponds to the observation in gel elec- 
trophoresis at high electric field (see Appendix lAl ITj|) 10] . 
Equation i|12|) provides a simple scaling estimate for the 
critical field of the abrupt transition on flow behaviours. 
For quantitative analyses, however, more details regard- 
ing the prefactor are required, which shall be discussed 
below. 

We consider cases with the electric field slightly lower 
than the critical injection field given by eq. 1)12(1 : Ei n = 
a E^ in 3) ( a < 1). Then there still remains a small free en- 
ergy barrier and an activation process of the permeation 
will be governed by the following free energy, which is 
a modified version of eq. 10 with the effect of outside 
events. 



F(y) - J' (AIM 2 - qeE in ^ dy 



(13) 



where All is the difference of osmotic pressures inside 
and outside of the slit and is given by 

An ~ k B T ( Jj - i) (w) 

Free energy eq. (|13|) takes a maximum value at 

y^D^l (15) 
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and the corresponding free energy barrier is 



AF 



3\2 



(16) 



Note that this description is valid at Ei n > E^ de ^ . At 
the threshold, F m = E^ de f \ the blob size becomes equal 
to the unperturbed coil £ ~ R\\ ~ R. At F in < E {de ^ 
the inside approach is valid. 

An important consequence of the chain deformation 
for E in > F< de ^ is that there are N/g blobs stuck at the 
slit entry. Every blob has an equal chance to climb the 
barrier toward a nucleation. The transition rate for this 
activation process is, therefore, estimated as 



__ X N 
1 9 



(17) 



where T\ ~ fi_D/£ is the mean escape time from the 
free energy trap, which is initiated from any single blob; 
fx ~ exp(AF/fc B T) with AF given by eq. {Tj|}. The 
factor -D/£ in t\ reflects the fact that the larger blob has 
larger contact area with the slit entrance^- The escape 
time is calculated by analysing the stochastic process of 
the barrier crossing. From eq. I|13fl . the Fokker-Plank 
equation for this process is 



dP(y, t) d 1 



dt 



dy i(y) 



dF(y) | d 
dy dy 



P(V, t) (18) 



where, 7(2/) ~ r\y is y-dependent friction coefficient. 
Evaluation for f\ is given as a mean first passage time 
from y ~ to y » y*[23j]. By assuming the bound- 
aries at y = and at y — yo >> y* to be reflective and 
absorbing, respectively, we obtain 



Wo 



dy y exp 



\k B T J 



v exp 
dy' 



f F(y')\ 

\ k B T / 



y'V(y') 



(19) 



where T>(y) ~ ksT/^{y) is the (^-dependent) diffusion 
coefficient of the partly confined chain|2J|. Following the 
standard procedures (utilizing the fact that the integrand 
of the first integral in eq. I|19|) is sharply peaked at y = 
y* , while that of the latter is very small around there), 
it is approximated as 



T\ ~ tjja 2 exp 



AF 
~k~RT 



(20) 



where to = rjD 3 / (ksT) is the characteristic time scale 
corresponding to the length scale of D. Substituting eq. 
(HHJl, (EDI) into eq. l(T7|). we obtain 



k = t D l w^ 3) (a, D)Nexp ( - 



AF 
~k~RT 



(21) 



with 



w 



(3), 



a, D) 



(- 



a e 



(22) 



At electric field below E^ de ^\ the chain deformation 
becomes irrelevant. The free energy barrier is, then, ob- 
tained through the inside approach (eq. 0J and the tran- 
sition rate is given by k = r-j" 1 ~ {f\D / R)~ 



-l. 



^V 2) (a, F)iVtexp 



AF, 



k B T 

(£j(per) < Em < E (def)j 



(23) 



with 



w 



(2) 



3 

a- 2 



(24) 



When the electric field is further decreased, the critical 
nucleation size of the hernia, y* n (eq. ©) becomes longer 
and equal to the length of the completely confined chain 

L = Na 5 ^ 3 /D 2 / 3 at F = F> er ) 



E (p er ) 



k R T ( D 



qeD 



- N- 



/ f|) 



(25) 



At electric field below F = F> er ), the free energy barrier 
is no longer AFi„ , but given by F„ (L) , which approaches 
the confinement free energy AFo (eq. Q) in the limit 
of zero electric field. The effect of electric field is, then, 
very small, and the transition rate is given by 



f^W 1 ) (a, D, AO exp 



F m {L) 
' k n T 



(F < F( per )(P6) 



with 



w {1) (a, D, N) 
a l L 




To obtain eq. (|26|l . one has to take notice that the lo- 
cation of the free energy maximum is no longer given by 
V = Vin ( cc l- ©1)7 but given by y = L (the length of the 
completely confined chain). In contrast to other regimes 
for higher electric field, the transition rate becomes ex- 
ponentially small with increase in the chain length (Sec. 

US- 
Hereafter, we refer to the migration dynamics at F < 

E (per)^ E (per) < E < E (def) &nd E (def) < £ < E (inj) ag 

regime I, regime II and regime III, respectively. 



V. ELECTROPHORETIC MOBILITY 

After having clarified the physics involved in the ele- 
mentary step of the device, it is now possible to discuss 
the mobility of the chain going through an artificial chan- 
nel consisting of a succession of geometrical traps. The 
mobility of the chain could be expressed as Q 



JL 

Pel 



to 



to 



(28) 



6 



The time lapse for a chain to travel an elementary unit is 
t = (neiEin)^ 1 {Lb + L D D/D), which is a controllable 
quantity from the device design (see Fig. P). Using eq. 
I)2ip. (|23|) and H2bj) . the explicit formula for the mobility 
is obtained; 



1.0 



JL ~ 



to 



t + AN~ 



to 



t + BN- 



(£(<W) < E < E (inj)) (29) 



< E < E (def)j ^ 



to 



t Q + (t D /wW{N))N- 



exp 



F,„(q 5 / 3 A7.D 2 / 3 ) 



{31) 



fe B T 

(E < E {per) ) 



where, for simplicity, we show the explicit dependence on 
N only. 

Figure 01 shows dependences of mobility on applied 
electric field for various chain lengths with different de- 
vice parameters (The bottom approximately corresponds 
to the experimental condition in ref. ( 5])). With the in- 
crease on field strength, a crossover is observed from a 
stuck state to a freely migrating state. Two trends are 
perceivable here: (i) If the chain is longer than a cer- 
tain threshold N*, a sharp crossover occurs around the 
well-defined field strength E*, which is almost length in- 
dependent. One can recognize, however, a clear trend 
that the crossover takes place at slightly lower field and 
becomes sharper with the increase in the chain length, 
(ii) For sufficiently short (but still larger than D) chains: 
iV(£>) = (£>/ a )5/3 < iv < TV*, on the other hand, the 
crossover becomes less sharp and is realized at much 
weaker field. 

Therefore, the dependence of the crossover field on the 
chain length is non-monotonic. As is clear from our anal- 
ysis, the second point, i.e., larger crossover field for longer 
chains, is a consequence of the weakly biased entropic 
trapping process (regime I), and the observation of the 
opposite dependence signals a strongly biased dynamics 
(regimes II and III). While the effect of chain length is 
exponential in the regime I, it shows up through a prefac- 
tor in regimes II and III. The threshold length separating 
above two opposite tendencies is given by 



N* =N { P er \E = E*) = ( ^ 



o 



(32) 



where a* = E* / E^ and the definition of is 
introduced below. 

Of great practical interest is how the electrophoretic 
mobility depends on chain length at a fixed dc electric 
field. An appreciable dependence on the chain length is 
expected if one tunes the electric field in place (slightly 
smaller than E^ m ^\ see Fig. 0J. This is exemplified in 
Fig. El for two different situations. For comparison, the 
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FIG. 4: Electrophoretic mobility of a polyelectrolyte in a de- 
signed channel as a function of electric field with the following 
parameters: (Top) D = 8 a, to = 5x 10~ 9 (sec), to = 0.2 
(sec) and N = 200, N = 500, 2 x 10 3 and 10 4 , (Bottom) 
D = 2 a, t D = 10" 6 (sec), to = 5 (sec) and N = 25, N = 50, 
250 and 10 3 . Around the transition region of a shown here, 
the dynamical behaviours correspond to regime II and regime 
III for the chain with (Top) N = 2 x 10 3 , 500, (Bottom) 50 
and (Top) N = 10 4 , (Bottom) N = 10 3 , N = 250, respec- 
tively. For the chain with (Top) N = 200, (Bottom) TV = 25, 
the transition is realized in the regime I. 



results obtained from the local friction picture (Appendix 
|A"|) are also shown in Fig. |SJ Note that in these figures, 
the right graphs approximately correspond to the exper- 
imental condition in ref. (0). 

The electrophoretic behaviours can be classified into 
the following four regimes. (O)There is no effect of the 
geometrical trap for very short chain (N < = 
(D/a) 5 / 3 ), then, the mobility is the same as that in the 
free draining limit. (I)For short chains (N^ < N < 
N^), the process is described by weakly biased en- 
tropic trapping (regime I). The mobility decreases ex- 
ponentially with the chain length (eq. lED)- Here, 
A^ er ) ~ {k B T/qeD){D/a) 2 / 3 E~ 1 from eq. (II)For 
chains with moderate length (N^ < N < N (def >), the 
permeation is driven by the electric field, but the confor- 
mation of the stuck chain is still in equilibrium (regime 
II). The mobility increases with the chain length accord- 
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FIG. 5: Length dependence of the electrophoretic mobility of 
polyelectrolytes at a given dc electric field with the following 
parameters: (Left) D — 8 a, to — 5xl0 -9 (sec), to = 0.2 (sec) 
and a = 0.068, (Right) D = 2 a, to = 10" 6 (sec), to = 5 (sec) 
and a = 0.075. Corresponding characteristic chain lengths 
are (Left) TV (D) = 32, TV (per) ~ 470,TV (de/) ~ 2800, TV C ~ 
1200, (Right) TV (D) ~ 3, TV (per) ~ 40, N {def) ~ 240, TV C ~ 270. 




10.0 



FIG. 6: Length dependence of the electrophoretic mobility 
of polyelectrolytes at a given dc electric field calculated ac- 
cording to the local friction model. The parameters are the 
same as Fig. Corresponding characteristic chain lengths 
are (Left) TV (D) = 32, TV (per) = N {ief) ~ 470, TV C ~ 2600, 
(Right) TV (D) ~ 3, TV (per) = AT( de/ ) ~ 40, TV C ~ 170. Around 
AT = TV (de/) , the dynamical regime is shifted from the almost 
unbiased entropic trap to the biased forced penetration. 




FIG. 7: A dynamical diagram of the electrophoresis in 
(TV, Ei n ) plane (Logarithmic scale). Lines separating regimes 



E 



(per) 



E {in3) (D/Rf /3 for the border between (I) and 
(II) and E (def) ~ E {ini) {D/R) for the border between (II) 
and (III). The dashed curve between white and gray areas des- 
ignates a crossover from a stuck to a flowing state, which is ob- 
tained from nto = 1 (see eq. 12810 . The gray area designates 
the regimes, in which the effect of the geometrical constric- 
tion is not apparent because of either too short chain length 
(left part) or sufficient strong electric field (upper part). At 
N >> TV*, this dashed curve is given by E* ~ [In TV] - (see 
eq. 1331 ') and this crossover exhibits a feature of a nonequi- 
librium phase transition in the limit of infinite chain length. 



tional behaviours of the stuck chain, which separates the 
regime II and III. We expect that a sufficiently long chain 
is markedly deformed at the entrance, which enhances the 
rate of the permeation. 

On the other hand, if one neglects the effect of 
hydrodynamic interactions (i.e., local friction picture), 
N ( per ) = N (def) is expectec } ; the regime (II) does not 
exist. Therefore, the deformation of the stuck chain is 
expected to be evident even for chains with moderate 
length. The computer simulation studies so far should 
be compared with this model[l5l liH ITtI ITs[ . 



ing to eq. JSH- Here, TV^-tt ~ (k B T/qeDE) 5 / 3 from 
eq. (0. (III)For longer chains (TV > TV (de/) ), the chain 
distortion at the entrance becomes apparent (regime III) , 
which leads to more enhanced mobility for longer chains 
(eq. |23)). 

As a brief summary, a diagram of the electrophoretic 
dynamics is shown in Fig. [7] It is worth mentioning a 
similarity between this diagram and the experimentally 
observed dynamical diagram of the gel electrophoresis 
[lOj (see also Appendix El EJ). 

The chain with TV > N^ per ^ is in the driven permeation 
regime, and the longer chains migrate faster. This is a 
consequence of the increased attempt frequency to climb 
the free energy barrier. How this attempt frequency de- 
pends on the chain length is a matter of the conforma- 



VI. DISCUSSION 

A. Comparison with Experiments 

First, let us compare our results with the experiment. 
In ref . [f| , electrophoresis of DN A was carried out using 
the narrow slit D = 90 nm. Then the sharp crossover 
from the trapped state to the freely migrating state was 
found around the electric field E* ve ~ 30 (V/cm), where 
E ave is the average electric field over the channel. From 
the device parameters used in the experiment, the corre- 
sponding electric field in the narrow slit is found to be 
E* ~ 50 (V/cm). The sharpness of the transition (de- 
fined as the ratio of the range of the electric field AE 
between zero and the asymptotic mobility to E*) is also 
found to be AE/E* ~ 0.6. Putting values correspond- 
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ing to DNA (segment length a ~ 60 nm and effective 
charge per segment q ~ 85) into eq. |Q , we obtain that 
£< in ^ ~ 20 (V/cm), and A£/£* ~ 0.5 is obtained from 
Fig. 0] (Bottom). Thus, we find a good agreement for the 
sharpness of the transition. On the other hand, we some- 
what underestimate the threshold field, since the experi- 
mentally observed value is slightly smaller than by 
a factor of a* ~ 0.1 (see Fig. 0J: Note that the precise 
value of a* depends on to, thus, the device parameters 
such as D, Ljy and L^. However, given the overall ac- 
curacy of scaling laws and the assumption involved in 
the simplified model (no leakage of the electric flux at 
the dielectric boundary, neglect of the effect of electro- 
osmotic flow), this discrepancy in the threshold field is 
not surprising. 



B. Long Chain Limit: Nonequilibrium Phase 
Transition 



It may be intriguing to examine dynamic behaviours 
at the asymptotic limit of the long chain length. As we 
shall argue now, a sharp crossover of the electrophoretic 
flow exhibits features of "phase transition" . We consider 
the situation, where the length of the device unit Lp, 
is large enough compared to the chain length. (Although 
this is, of course, not the case for very long chains, one 
can take the limit, in which the device size is made ar- 
bitrary large in proportion to the chain length.) In this 
limit, the threshold field separating a stuck and a flowing 
state becomes 



E* 



[In AT] 



and the electrophoretic mobility is expressed as 



1 



Me; 



1 + const. N e 



(33) 



(34) 



where e ~ (a* — a) /a is a reduced electric field. The 
crossover becomes shaper and sharper with the increase 
in N and becomes singular in the limit of infinite chain 
length as announced. Therefore, this dynamical state 
change (from the state with zero to some finite current) 
caused by the external field is formally regarded as a 
nonequilibrium phase transition. 

As another possible choice of the limit, one may con- 
sider increasing the chain length for a fixed device size. 
In such el CclSG, <X chain interacts with many elementary 
device units simultaneously, and some unscripted new 
features should appear (l8| . 



C. Practical Aspects 

A unique property of the present system is that there 
is a regime where longer chains migrate faster. This is 
a consequence of the larger contact area with the slit 



entrance for larger chains For long enough chains 
(JV(<fe/) < TV), or under the strong electric field (E^ def ^ < 
E), this effect is further enhanced by the electric field 
induced chain distortion at the geometrical trap, which 
generates more blobs for longer chains as steady state 
conformations (regime III). An increase in blob numbers 
increases the frequency for the trial to climb the free en- 
ergy barrier as is clearly seen in eq. (|21|l . If one neglects 
the effect caused by the solvent flow, the regime II dis- 
appears. Then, the chain stuck at the slit entrance and 
driven by the electric field is inevitably distorted. This 
explains the reason why the stuck chain deformation was 
observed in Monte Carlo simulations^^, while it is not 
apparent in real experiments using moderate length of 
DNA0. 

For very long chains of N >> N c — A /to (or N c = 
(B/t ) 5/3 if (B/t ) 5 / 3 < N^), however, the chain 
trapping becomes irrelevant (to >> k _1 ), and the mo- 
bility approaches to the free draining limit. This bounds 
the range for efficient separation in a single run of elec- 
trophoresis at a given set of parameters. However, we 
now know the trend on how this appropriate length range 
depends on controllable parameters such as size of chan- 
nels (D, Lr>, Lp) and electric field. The same also applies 
for the migrating speed of chains, which is another impor- 
tant factor for the practical purpose. In addition, such 
microfabricated devices could be easily assembled into 
practical integrative systems. Therefore, by construct- 
ing multiple channels with designed different conditions, 
which operate in succession or in parallel, it may be pos- 
sible to achieve an efficient separation over a wide range 
of chain length. This is, in fact, of advantage to use de- 
signed channels for separations. 



VII. FINAL REMARKS 

Our present analysis is based on the steady state rhe- 
ological behaviours of the chain. This approach is moti- 
vated by the microscopic images of the experiment[^,|i|, 
in which the transient behaviours (the relaxation after 
reentering the thick region, the stretching under high 
electric field gradient, etc.) seem not to be important 
very well. However, if the period of the device is made 
shorter, successive injection and ejection events start to 
correlate, and the transient dynamics of the chain be- 
comes relevant [l5l Il7l Il8j . Also, recent experimental and 
simulation studies have observed the transition between 
two distinct states with different migration speeds using 
relatively deep channels, in which DNAs may be trapped 
deep in the well and can not escape from it [17], [l|| ■ Such 
features are not studied here, but should be incorporated 
in the future analysis. 

In summary, we have shown that the deformation of 
the stuck chain and the effect of the flow are important 
ingredients for the DNA electrophoresis in the designed 
channel under consideration. Such rheological proper- 
ties in a single chain level are expected to manifest their 
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importance in many other situations in the dynamics of 
polymers in restrictive environments. 
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independent electrophoretic mobility. This corresponds 
to the observation in gel electrophoresis that the entropic 
trapping effect disappears at electric field higher than a 
certain threshold value [lOj. (Note that the topological re- 
striction imposed by gels with rather large pore size may 
be better modelled by the cylindrical pore geometry, for 
which we find the same scaling form of the injection field, 
see Appendix IbI) 

The activation process at a{= Ei n / Ey^j') < 1 can be 
analyzed similarly. Now the free energy (eq. i|13|) ) f° r 
the permeation process is modified as 



APPENDIX A: ELECTROPHORESIS BASED ON 
THE LOCAL FRICTION ASSUMPTION 

In this appendix, we briefly discuss the electrophoretic 
dynamics based on the local friction picture, where the 
role of the solvent flow is neglected. Although this does 
not correspond to the real experiments, it provides useful 
information to analyze the result from computer simula- 
tions performed so far. In addition, this argument may 
be relevant to the gel electrophoresis, where the hydrody- 
namics interactions are screened at the scale of the pore 
size (see also Appendix |B|) . 

Then, the rheology of the stuck chain is modified as 
follows. Instead of eq. (jSJl, the electric force acting on 
the stuck chain is 



f(R) 



geE m DN = geE in DR* 

R 0,3 



(Al) 



Here and in the followings, the subscript (R) indicates 
the local friction (Rouse) dynamics. One finds for the 
deformation field 



qeD 



(A2) 



At electric field above E^J' , the stuck chain is deformed, 
and its conformation is a parallel alignment of blob of size 
(cf. Fig. EJ. 



The length of the stuck chain along the entrance if 

2 

/ qeE in D\ 5 



R. 



(R)\\ 



aN 



V k B T 



(A3) 



(A4) 



The condition for the smooth chain injection is D > £(_r), 
which leads to the critical injection filed; 



E, 



(inj) 
(R) 




(A5) 



At Em > , the trapping at the slit entrance is irrel- 

evant, which results in the length and also electric field 



F(R)(y) 



(An 



(R) 



D z -qeE ln g D ^)dy (A6) 



with 



This produces a free energy barrier 



AF, 



(R) 



(1 -a 9 / 5 ) 2 



k B T 



and its location 



D 



{!-<**) 



(A7) 



(A8) 



(A9) 



By taking into account that the diffusion coefficient of 
the partly confined chain is now given by T>r R \(y) ~ 
(kBT/rj)(a/D) 2 / 3 y~ 1 , the transition rate is calculated as 



4.-1 (3)/ n\Ar ( ^ F (R) 

K (R) = t D w ( R )( a ' D)Nexp I — -j—^- 

F (def) F E (inj) 



(A10) 



with 



(All) 



Note that this description is valid at Ei n > E^^ 

(regime III). At the threshold, E in = e9^/\ or 

equivalently a ~ (D/a) 5 / 3 N~ 1 , we have yf R -, ~ 

(a 5 / 3 /D 2 / 3 )N — L, which is the length of the com- 
pletely confined chain. Therefore, at this point, we 
have a crossover to the almost unbiased entropic bar- 
rier regime (regime I) and the regime II does not exist 

{E^^ ~ EygP)- The transition rate in regime I is cal- 
culated as 

K = rtwV(a, D, AOA^expf-^^V^) 



k B T 

Ei„, < E, 



(def) 
(R) 
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where F in ^ describes the free energy of the permeation 
process based on the inside approach, and 



(1) 

(R) 



(a, D, N) 



5/3 



L 



D-aL 



D 



D-aL 



1 

a 



The dependence of the electrophoretic mobility on the 
chain length is represented in Fig. H3 



where R^ de ^> is the overall size of the deformed stuck 
chain and 7 fl (d C /> is the friction coefficient of the outmost 
blob. 



f kgTa 
\qeE m D 2 



R (def) 



V Bin 



i?( de {B4) 



The force acting on the inner blobs (r < i?( de ^)) m a 
radial direction is the sum of drag forces acting on outer 
blobs; 



APPENDIX B: CHANNELS WITH 
CYLINDRICAL PORES 

If the channel is composed of the alternating wide and 
narrow cylindrical pores, instead of deep and shallow 
slits, the electrophoretic dynamics of the chain is altered. 
In this geometry, the convergent field near the entrance 
of narrow constriction is E(r) ~ Ei n D 2 /r 2 . The stall 
force for the unperturbed stuck chain is evaluated as 



/ ~ 'foUelE = 



qeE m D 2 
aR 



(Bl) 



where 70 — rjR is the friction coefficient for an unper- 
turbed chain, and E ~ E(r)r 2 / R 3 dr is the average 
electric field acting on the stuck chain. By comparing 
this force with fc^T/i?, the deformation field is obtained 



E (def) 



ksTa 
qeD 2 



0). 



(B2) 

E(def) does 



In contrast to the slit geometry (eq. 
not depend on the chain length. 

At the electric field above E^ de ^\ the chain is highly 
compressed toward the narrow pore and the resultant 
conformation can be viewed as a dense stacking of blobs, 
whose size shrinks with approaching to the pore entrance 
i.e., an inhomogeneous semidilute solution (Fig. [SJ). The 




FIG. 8: A stuck chain at the entrance of the cylindrical pore. 
The colour strength represents the segment density. 



size of outmost blobs is obtained from 

k B T 



lR{d . f) H el E(R(^f)) 



(B3) 



f(r) 



qt E in D 2 (\ 1 

r i?( de /) 



dr 



(B5) 
(B6) 



By balancing this force with fcgT/^r), the blob size at 
r > D is expressed as 



k B Ta rR( de ^ 
qeE in D 2 R( de f) - \ 



(B7) 



At r < D, the blob does not further shrink any more; 
£(r) = £(-D) at r < D. The size of the deformed chain 
Rl de f) i s determined from the segment number conserva- 
tion; 



9(r) ,, 



d 6 r ~ N 



(B8) 



where gir) = (£( r )/ a ) 5 ^ 3 is the number of segments in a 
blob at r. 

Since D « R( de $\ the blob size at the entrance (r = 
D) is readily approximated as 



ksTa 
qeE in D 



(B9) 



We expect the smooth injection when £(D) < D. There- 
fore, we find that the injection field £K lT w is sparsely 
larger than the deformation field, but no dramatic dif- 
ference is expected; E^ m ^ ~ E^ de ^ . Thus, the injection 
field is almost the same as that in a slit geometry (eq. 
iQlJl). It means that once the stuck chain is deformed at 
the entrance, there is almost no free energy barrier left, 
and the smooth injection takes place (the range for the 
regime III is very narrow), therefore, the use of the in- 
side approach (Sec. IIV Afl is justified in this case. The 
field E^ er \ below which the process is described by the 
weakly biased entropic trapping, is the same as that in 
the slit geometry (eq. H25(l h 

The case with local friction picture can be worked out 
similarly. For a chain obeying local friction (Rouse) dy- 
namics, we find that 



E, 



(def) 
(R) 



ksTa 
qeD 2 



N~ 



(BIO) 
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Close to the entrance (r « R^ de ^) 
approximated as 



the blob size is 



( k B T 

\qe,E in D 



(Bll) 



From this, we find that the injection field is given by the 
same scaling form as that in a slit geometry (eq. (1A5|) S ) 

The field , below which the process is described by 

the weakly biased entropic trapping, is the same as that 



p (per) 



E 



(de/) 
(R) 



in the slit geometry (eq. (|A2j) : note that 
for a slit geometry). 

In contrast to the slit geometry, the electrophoretic 
mobility does not depend on the chain length in the 
driven biased activation regime (regime II and III) irre- 
spective of the presence or the absence of the flow effect. 
The length dependence is expected to show up in the 
regime I only, where the longer chains migrate naturally 
much more slowly. In gel electrophoresis, this entropic 
trapping process is observed, when the pore size is rela- 
tively large and electric field is rather weak. Here, one 
comment might be useful. In gel electrophoresis, the en- 
tropic trapping is replaced by the reptation process at 
the electric field higher than a certain threshold. Experi- 
ments show that this threshold is lower for lower gel con- 
centrations, but seems not depend on the chain length 
of the migrating polymer, which accords with the fea- 
ture in our injection field 10]. Of course, this comparison 
should be regarded as very qualitative, but the trend is 
apparently the same. 



we breifly note the connection between electric and flow 
induced injection problems. The control parameter in the 
latter case is the current of the solvent flow J (m 3 /sec) 
through the pore. Near the entrance of the pore, the 
flow is convergent, and the flow field is calculated as in 
the same way as the electric filed in Sec. lIVBI fslit]) and in 
Appendix (capillary). The characteristic currents can 
be derived using the hydrodynaic-electric equivalence. 

In the slit geometry, we measure J as the solvent cur- 
rent through the slit per unit length. From eq. I|25|l . Q 
and l|12fl . we find respectively, 



j(per) 



ksTfD^ N _, 
r/a \ a 



j(def) 



j(inj) 



k B T 



-N- 



k B T 
■qD 



(CI) 



(C2) 



(C3) 



In the capillary geometry, we measure J as the overall 
current through cross sectional area of the pore. From 
eq. (23 and l|B2"|) . we find 



j(per) 



k B T ( D 



N~ 



(C4) 



APPENDIX C: CONNECTION TO THE 
FLOW-INJECTION PROBLEM 

Polymers can be driven into a narrow slit or a capillary 
by applying a flow of the solvent pjj lU |22, |2g. Here 



j(def) „ jiinj) „ k JL (C5) 
V 

The critical injection current in eq. (|C5|) is consistent 
with the result obtained by different approaches |2ll l2q . 
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